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Abstract. - We discuss the quasiclassical approximation for the equations of motions of a 

O^ ' nonlinear chain of phonons and electrons having phonon mediated hopping. Describing the 

phonons and electrons as even and odd grassmannian functions and using the continuum limit 

we show that the equations of motions lead to a Zakharov-like system for bosonic and fermionic 

^^ ' fields. Localised and nonlocalised solutions are discussed using the Hirota bilinear formalism. 

Nonlocalised solutions turn out to appear naturally for any choice of wave parameters. The 

bosonic localised solution has a fermionic dressing while the fermionic one is an oscillatory 

localised field. They appear only if some constraints on the dispersion are imposed. In this 

case the density of fermions is a strongly localised travelling wave. Also it is shown that in the 

multiple scales approach the emergent equation is linear. Only for the resonant case we get a 

nonlinear fermionic Yajima-Oikawa system. Physical implications are discussed. 
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^^ . Introduction. - The idea of treating classically systems with fermions goes back to sev- 

'^ ' enties to the papers of Berezin [1] , Casalbuoni [2] . This is done using grassmannian functions 

(elements of a Grassmann algebra). A simple representation can be done using, for instance, 

combinations of Dirac matrices, namely a = 71 + 172, /3 = 73 + i74. The algebra generated 

by (1, a, /?, Qf/S) is a grassmann algebra with two generators. The even sector (bosonic or 

S ' commutative) is given by linear combinations of (1,q;/3), while the odd sector (fermionic or 

^ , anticommutative) is given by linear combinations of (a, /3). 

\rl \ Today, grassmannian description is used intensively in mathematical physics. It is a very 

'Oj . interesting idea because it could be implemented to nonlinear evolution equations containing 

C^ ' bosons and fermions, for which the quantum description is difficult to do. Moreover, since 

the quantum mechanics deals mainly with observables and spectra, classical description offers 

an idea of the dynamics of solutions for the equations of motions. In this view there are 

some results related to the dynamics of solitons in integrable equations containing bosons 

and fermions but linked supersymmetrically [3]. With the aid of a supersymmetric bilinear 

formalism these type of equations have been solved. Interesting enough the existence of 

fermions makes the soliton interaction to be nonelastic having a dressing in the fermionic 

part. 
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Due to these new results it seems naturally to try to extend them to some physical models 
involving bosons and fermions and to sec how the dynamics of the localised solutions is 
influenced by the fermions. 

The main realm is of course the problem of localised structures in biomoleculcs (proteins, 
DNA, etc.). Since the seminal ideas of A. S. Davydov [4], it has been continuing an intense 
quest for finding experimentally these nonlinear structures (solitons, breathers, etc.) [5]. In 
any case, insofar the results are controversial. 

In this paper we consider the simplest model namely a nonlinear chain and describe the 
phonons and electrons as functions having values in the even and the odd sector of an infinitely 
generated grassmann algebra (this can be done either directly or using coherent states). The 
choice of SSH hamiltonian is not only for simplicity, but also it captures the essential informa- 
tions when one deals with a nonlinear lattice having a coupling between charge and structure 
(which explains charge conduction, polaron formation etc.) [6], [7]. 

In the continuum limit the differential-difference equations of motions become a system 
which looks like the classical Zakharov system [8] with fermions. Interesting enough this 
system admits exact bilinearisation in the Hirota formalism and two type of solutions. One 
type is exponential and nonlocalised although it appears to look like a kink type solution. But 
the presence of fermionic parameters at denominators breaks the localisation. The other type 
is localised even though is dressed with fermionic parameters as well. This type of solutions 
is responsible for a soliton-like propagating structure of the fermionic density in the chain. 

Equations of motion. - The hamiltonian we arc working with is 

H = ^ -^ + W{Un+l - Un) + (cl+iCn + c]^C„+i) [-tg + 7("n+l - "«)] (1) 

n 

where Un is the displacement from the equilibrium position of the particle n and c|j(c„) are the 
creation (annihilation) operators for fermions on the site n. Also the potential W(un+i — Un) 
can be harmonic as in the Su-Schriffer-Heeger (SSH) problem [6] or anharmonic. We do not 
consider the spin indices because we are not interested here in magnetic phenomena. The last 
term which contains the phonon mediated hopping is responsible for the interesting physics of 
the system. The classical hamiltonian can be obtained by averaging the quantum one using 
the product of the following coherent states (or directly by substitution) 

|x„>=e-X"X"/2e'^"^"|0>, 

where Xn{t) are grassmann- valued odd (anticommuting) complex functions. Also M„(t) are, 
from now on grassmann-valued even (commuting) functions (they are not real functions any- 
more) . They belong to the even and odd sector of an infinitely generated grassmann algebra. 
Complex conjugation is defined such that (Z1Z2)* — Z2z\ in accordance with the conventions 
in [9]. As we are dealing with the classical case we assume further that all fields and their 
derivatives commute or anticommute, depending in the usual way on the bosonic or fermionic 
nature of the fields. 

First we consider the harmonic approximation, namely Wiun+i — Un) = K /2(un+i — Un)"^. 
After averaging on the coherent states we find the following classical hamiltonian. 

^ = X! Y" + y (""+1 ^ "")^ + (Xr*+lX« + XnXn+l) [-to + l{Un+l " ««)] (2) 

n 

The equations of motions are given by [2] : 

M„ = K{Un+l + Un-1 - 2w„) + ■j{Xn+lXn + XnXn+1 - X;^X«-1 ^ Xn-lXn) (3) 
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iXn = -toiXn+l + Xn-l) + 7 [(^n " "n-OXn-l + ("n+l - '"n)Xn+l] (4) 

A tractable version of the above equations can be obtained using their continuum hmit. 
Considering x — en we keep only terms up to order two in the expansions, i.e. w„±i = 
u(a;, i)±M;j^+Uj;2;/2 and the same for Xnii- After gauging the fermionic field x — > xexp(— 2itot) 
and calling K = c^ we find the following system: 

utt - c^u^x = 27(x*x)a; (5) 

iXt + toXxx - '2'lUxX = (6) 

As expected we get the Zakharov system with bosonic and fermionic components. The corre- 
sponding hamiltonian for these equations is given by: 

H = j dx{-ul + —u\ + 27x*XMx + t^xlXx j (7) 

Localised and nonlocalised solutions. - In order to solve the equations of motions |(SJ), 
© one can write all the quantities in terms of the generators of Grassmann algebra and then 
solve the equations layer by layer [10]. Even though the results can be very interesting, they 
depend on the number of generators. A different approach which is more appropriate for 
computing the soliton-likc solutions is the Hirota bilinear method [11]. The idea behind it is 
to implement a nonlinear substitution for the field such that the nonlinear system is turned 
into a bilinear one having the form of the dispersion relation when is written using the Hirota 
operators. On the other hand the nonlinear substitution somehow swallows the singularities 
of the solutions of nonlinear system and the bilinear one become a system which contains 
holomorphic functions. In the case of integrable systems this fact is rigurously true. But 
in the case of nonintegrable systems this works only for special solutions which, because of 
the holomorphic character, can be written in terms of exponentials. This strongly alleviates 
all the calculations. On the other hand the adaptation of Hirota formalism to grassmannian 
functions is straightforward. 

The nonlinear substitutions are: 



X{x,t)^—-^ — , u{x,t) = — -d^logF 

2'y t 7 

where G{x, t) is an grassmann anticommuting complex function and _F is a grassmann com- 
muting real function. Plugging into the nonlinear equations we find the following bilinear 
system: 

{Til - c^T)l)F • F = GG* (8) 

(zDt+toD2)G.F==0 (9) 

where D"/ • g = [d^Y f[x + e)g{x — e)\e=o is the Hirota bilinear operator. Now, usually the 
second bilinear equation gives the dispersion relation and the other one is only a constraint 
on the phases. The 1-soliton solution usually is given by the following ansatz, 

G^Ce", ^=1 + 6"+"'+^ (10) 

where (^ is a complex anticommuting grassmann parameter (^ = C-R + ^C/)j V — kx+ujt where k 
and u) are commuting complex grassmann parameters and a constant phase. After plugging 
the ansatz into the equation we find the following relations for the parameters: 

[{u; + u;r-c'{k + kr]e^ = lcC (11) 
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iu + tak'^ = aCC* (12) 

The second equation and its complex conjugate arise as overall factors multiplicating terms 
containing ( and (* . Accordingly these factors must be zero or proportional with the product 
(^^* by a complex factor a. 

Case a = 
As we pointed out, the eq. (|12|l gives the dispersion relation of the soliton and the equation 
Hll|l is necessary to compute the phase cj). Indeed, in our case, 



2(w + tJ*)2-2c2(fc + fc*)2 

i.e. is proportionally with the product of zetas. So, there are absolutely no restrictions on k 
and (. For simplicity let's call the denominator denominator A~^. Our 1-soliton solution will 
be (up to the factors) 

G Ce" , ,, F^ (fc + fc*)ACCe''+''* 

^v ' ^ p l + A((*eV+n" ^ ' ' F 1 + ACC*e''+''* 

In the grassmann algebra we can compute exactly the inverse of an commuting object, namely 
l/(a + 6C1C2) = 1/a — C,iC,2b/a? and this is because of the fact that Q — Q. Applying this rule 
to our soliton we shall find that our soliton is not at all a soliton but a pure nonlocalised 
exponential solution. 

X{x, t) = Ce", u{x, t) = {k + fc*)ACC*e''+''' (13) 

One can see that the bosonic solution is purely grassmannian (its square is zero), so we can 
speculate that it exists only at the quantum level. Accordingly, the usual method of finding 
localised solutions for Hirota bilinear equations here fails due to the grassmannian character 
of the parameters. 

There is a big drawback of these solutions. They have no physical meaning since they 
grow exponentially and the energy as well. 

So we are forced to look for a method to find localised solutions. A posibility is to treat 
diferently equation Hll|l . namely to take it as a dispersion relation while the equation H12|l 
will give constraints on the solitonic parameters. For instance we can consider that e"^ is not 
at all proportional to the product of zetas but just a pure number (let's say 1/8) and zeta's 
enter in a the relation of w and fc i.e 

(w + uj*f - c2(fc + fc*)2 = e-'^CC72 

If we put lli = uijf + iojj, k ^ kjf + ikj and e*^ = 1/8 we find the following relation 



^- - \^^^^+CC = cfc. + ^CC* = ckn + -L-OCn 

the last equations being consequences of the grassmannian character of ( and (*. Also the 
last equation is a real grassmannian quantity despite the i factor because of the definition of 
complex conjugation (^1^2)* = 2^2-^1. Plugging this relation into the equation (|12|) one obtains 

CC \ .^2^2 CC 



'^-"'.^ + i^^j' "^-'"r^"i^ MF 



O'^R 
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These are the constraints in the sense that kj and w/ are no longer free. So, instead of having 
four free parameters kji, kj, C^, Q we have only three fc/j, Qn and C/- With all these relations 
one can write (after some algebra) completely the 1-soliton solution, 

u{x,t) = '^sech% fe^* + T^e*«*) (14) 

2 V KflC / 



X{x,t) = -sechfjoexp 



A—-X - i{k%to - -—)t 



(15) 



where rjQ = k^x + ckpt. One can see that the fermionic field x is indeed a localised field 
with an oscillating amplitude. The problem of localisability for the bosonic field is not so 
obvious because of the fermionic dressing in the last exponential. In any case expanding the 
exponential (taking into account that the square of C,C,* is zero) we see that the problematic 
part of u{x,t) is {t/kjicjsech'^rjQ which has finite limits as time goes to infinity. So we have 
indeed a localised solution which is expressed using bosonic and fermionic parameters. The 
bosonic field contain a fermionic correction and the fermionic field is a soliton with oscillating 
complex amplitude. Here the physical interpretation is clear. This can be seen as a classical 
analog of the localised polaron wave function for the SSH problem [7]. With our method we 
can in principle prove the existence of localised states. Using this solution one can compute 
the behaviour of the fermionic density, namely cJjC„. In the continuum limit this expression 
is given by 

X*{x,t)x{x,t) = —sech?-qo 

which shows the strongly localised character of the charge density travelling wave. Plugging 
these solutions in the expression of the hamiltonian we find the energy: 

One can see that the dressing of the bosonic solution and the fermionic field gives a purely 
grassmannian correction. 

Case a ^ 
In this case, considering a — j-i + iv the relations for fc/ and uji have the same structure. Only 
the factors of zetas are translated with some terms proportional with /i and u. 

Of course, a big advantage of this model is that the underlying bosonic part is just a linear 
Klein-Gordon equation, so we do not have solitonic phenomenology. We can do absolutely 
the same machinery for a nonlinear completely integrable bosonic equation with rich solitonic 
phenomenology. For instance, we consider a fermionic extension of the Boussinesq equation 

Utt - Uxx - QUxUxx - Uxxxx = {XX*)x (16) 

iXt + Xxx - 2uxX = (17) 

This system can be seen as a continuum version of the equations of motions for the hamiltonian 
(2.1) with the potential (up to coefficients) VF(u„+i— u„) = l/2(u„+i — w„)^ + l/3(u„+i — u„)^ 
for a weak phonon mediated hopping (we neglected the coefficients in the equations). Then 
with the same nonlinear substitution we shall get: 

(D2-D2-D^)F.F = G*G, (*D, + D2)G.F = (18) 
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The ansatz for l-soliton solution gives the foUowing dispersion relations 

[{uj + uj*f-{k + k*)^~{k + k*)^]e^ = ^C<:, iiu + k^=aCC (19) 

One can see that everything works in the same way except that A^^ is translated with (k+k*)'^. 

Multiple scales approach. - Even though the continuum limit used above is a naive one 
we can prove that using the systematic multiple scales method [12] one gets the unidirectional 
approximation of the system (2.5), (2.6). 

In order to do this we assume that the competition between dispersion and nonlinearity 
occurs at large scales of space and time. Accordingly, we introduce the following stretched 
variables: 

f = e(n — Vgt), T — e^t 

and consider that the fermionic complex field is just a slowly modulated fermionic wave. 

In the exponential, w is exactly the dispersion relation of the second discrete equation from 
the system (3), namely lu — — 2io cosfc. Also Vg in the definition of the streched space is the 
group velocity corresponding to this dispersion. For the bosonic field we take, 

a.„(i)=5]e^t¥,(e,r) 
i>i 

Using the discrete equations 10} we shall find at the order e^ 

j$iT- = —uj2^i^^ — 27Cosfciyi^$i 
where W2 = dPu/dk'^. Also from the first discrete equation Q we shall find 

e3(i;2 - c^)Wi^^ + 0{e^) = 276^ cos fc($*$i)5 + 0{e^) (20) 

where c^ is the phononic velocity. Assuming that we have the nonresonant case namely Vg ^ c 

we find 

27C0sfe^^^ 
^1? = 2 T*t*i 

Vg — C^ 

Introducing this in the above equation and taking into account that the square of the $1 is 
zero we get a linear equation 

i<^ir + t^2$i« = (21) 

Accordingly, in the nonresonant case there is no localised solution. 

For the resonant case i.e. Vg — c we have to change the scaling of the fermion field 

X„(t) = e^('="-"*)^e^+i/2<3>,(C,r) 

which leaves unchanged the second equation but the first one will be 

e\vl - c2)iy2« - 2e%iyi5, + ^(e') = '2e^lcosk{<^l<^^)^ + 0{e^) 
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Because Vg — c, finally we have the following system: 

Wir + — cos /c($i*i) = (22) 

-"a 

i$ir + W2$i^5 + 27COS kWi^^i = (23) 

which is indeed the uni-directional approximation of the Zakharov-type system, the so called 
Yajima-Oikawa system [13]. This equation has also a bilinear form and we can compute the 
solutions in exactly the same way. 

Conclusions. - The dynamics of localised solutions in a nonlinear chain with fermions 
described by SSH model, is solved here using bilinear formalism adapted to a grassmannian 
formulation of the equations of motions. The main problem here was that the usual way of 
finding solitonic solutions gives purely exponential nonlocalised solutions. This fact happens 
even in the case of supersymmetric integrable equations [14] using the celebrated Darboux 
transform. Anyway including the grassmannian parameters into the dispersion relation we 
have found really localised solutions but with some constraints on the parameters. So, one 
can speculate that in the nonlinear equations with fermions localised solutions exist only in 
special cases. 

Also from the above multiple scale approach we can see that the existence of localised 
solutions is somehow problematic at least for SSH model inasmuch the general nonresonant 
case in the standard scaling goes to a linear equation for the fermionic field. Of course one can 
argue that a totally different scaling of space, time and fields might give a nonlinear equation 
for the fermionic field. We do not believe this since purely fermions cannot selforganize into 
solitons due to the exclusion principle. In order to appear localised solutions into a fermionic 
equation this one must be coupled with a bosonic one. 
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